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ABSTRACT 

The  stability  of  a  circular  cylindrical  shell  subjected 
to  a  moving  ring  load  with  a  constant  velocity  has 
been  examined  in  detail  when  both  longitudinal  and 
transverse  inertia  effects  are  included,  using  both 
the  Timoshenko  and  the  Fltigge  equations.  A  lower 
resonance  speed  has  been  found,  which  apparently 
has  been  concealed  in  a  recent  paper. 
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NOMENCLATURE 


a 


a,  b 


D.  A. 


E 

h 

P 

PO 

t 

T 


u,  w 
u,  w 
V 

VCR 

v 

VCR 

v0'  V1 
x 

X 

x0,  Xj 


mean  radius  of  shell 
see  Eq.  (27) 
see  Eq.  (28) 
see  Eqs.  (61)  and  (62) 
see  Eqs.  (39)  and  (40) 
see  Eq.  (34) 

Young's  modulus 
shell  thickness 
(1  -  *2)p/\E 

ring  load  per  unit  circumference 

at /a 

time 

u/a,  w/a 

axial  and  radial  displacements,  respectively 
velocity  of  ring  load 
physical  resonant  speeds 
V/a 

nondimensionalized  resonant  speeds 

see  Eq.  (20) 

x/a 

axial  coordinate 
see  Eq.  (52) 
x  -  vt 
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2 

al* 


2 

a2 


Ql'  a2 
32 


A 

A' 


6(z) 

c 

( 

\ 

l> 

« 

P 

<P 

n 


=  E/p(l  -  „2) 

=  see  Eqs.  (73)  and  (74) 

=  see  Eq.  (75) 

=  1  -v2 

=  see  Eq.  ( 1 3) 

=  see  Eq.  (49) 

=  Dirac  delta  function 

=  \hj3 

=  +VTT 

=  h/a 

=  Poisson's  ratio 
=  Fourier  transform  variable 
=  shell  density 

=  v/V7 

=  Fourier  transformed  quantity 
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SECTION  I 
INTRODUCTION 


In  recent  years,  the  problem  of  the  dynamic  response  of  circular 
cylindrical  shells  has  received  considerable  interest  because  of  their  use 
in  aerospace  vehicles.  In  particular,  many  investigations  have  examined 
shells  subjected  to  traveling  loads:  Nachbar  [  l]*  considered  the  dynamic 
response  of  an  infinitely  long  cylindrical  shell  to  moving  discontinuous  loads. 
Brogan  [2]  considered  the  problem  of  a  uniform  pressure  front  moving  over 
a  finite  cylindrical  shell  which  was  taken  to  be  internally  pressurized,  but 
the  resonant  modes  were  not  considered  and  the  solution  obtained  was  not 
valid  for  all  values  of  time.  Bhuta  [  3]  obtained  the  transient  response  of 
a  finite  shell  which  was  not  pressurized,  but  the  resonant  modes  were 
studied  in  detail.  It  was  also  pointed  out  in  [  3]  that  a  critical  velocity  exists 
for  a  semi -infinite  shell  for  which  the  deflections  become  unbounded.  The 
same  critical  velocity  was  obtained  by  Prieskin  [4]  for  an  infinitely  long 
cylindrical  shell. 

In  all  of  the  forementioned  investigations,  only  the  axisymmetric 
loading  was  considered.  However,  in  [  l],  viscous  damping  in  the  radial 
direction  was  assumed  to  be  present.  In  [2],  [3],  and  [4],  the  investiga¬ 
tions  concerned  the  response  of  undamped  shells.  Also,  [  l]  differs  from 
the  rest  in  that  the  acceleration  in  the  axial  direction  was  considered, 
whereas  in  other  prior  investigations  only  radial  inertia  was  taken  into 
account.  The  critical  velocity  of  [4]  due  to  a  bending  resonance  does  not 
appear  in  the  solution  of  [  1  ] ,  but  rather  a  much  greater  critical  velocity 
due  to  axial  motion  resonance  appears.  The  resonance  due  to  bending  seems 
to  have  been  concealed  because  of  the  assumed  presence  of  damping. 
Unfortunately,  however,  as  is  shown  in  the  present  investigation,  the 
resonance  due  to  bending  occurs  at  a  much  lower  velocity  than  does  the 


1 


Numbers  in  brackets  designate  References  at  end  of  paper. 
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resonance  due  to  axial  motion  and  is  important  for  aerospace  vehicle  design. 
The  resonant  character  of  the  solution  was  brought  out  also  in  [  3]  by  having 
let  the  shell  length  become  successively  larger  and  having  noted  that  at  the 
critical  velocity  the  dynamic  loading  factor  became  infinite. 


o 

SECTION  II 

STATEMENT  OF  PROBLEM  AND  METHODS  OF  SOLUTION 


It  is  the  purpose  of  this  paper  to  investigate  the  dynamic  stability  and 
solve  the  problem  of  the  dynamic  response  of  an  infinitely  long  cylindrical 
shell  subjected  to  a  ring  load  moving  with  constant  velocity,  taking  into 
account  the  accelerations  in  both  the  radial  and  the  axial  directions.  It  is 
shown  that  even  when  both  accelerations  are  considered,  the  resonance  due 
to  bending  occurs  at  a  lower  velocity  than  the  one  due  to  axial  motion. 

The  problem  is  solved  using  Timoshenko's  equations  [5],  as  in  [  l] 
but  without  the  presence  of  damping,  and  Flugge's  equations  [6].  The  results 
of  the  analyses  are  compared,  and  it  is  found  that  Timoshenko's  equations 
suffice  for  values  of  the  load  speed,  up  to  the  first  critical  speed.  It  is  found 
that  above  that  speed,  there  is  really  no  proper  or  physically  meaningful 
steady- state  solution,  as  was  concluded  also  in  [4],  and  that  an  initial  value 
problem  should  be  solved  for  such  speeds.  Near  the  axial  resonance,  the 
effects  of  shear  deformation  and  rotatory  inertia  should  be  included  also. 

In  addition,  expressions  for  dynamic  amplification  in  the  deflections  are 
derived. 

The  analysis  is  carried  out  using  Fourier  transforms.  The  resonances 
or  instabilities  are  shown  by  the  fact  that  for  certain  velocities  the  inversion 
integral  fails»to  exist. 
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SECTION  III 
ANALYSIS 


Solution  Using  Timoshenko  Equations 

The  shell--of  thickness  h  and  radius  a,  and  of  a  material  having  a 
modulus  of  elasticity  E,  Poisson's  ratio  v,  and  density  p--is  loaded  by  a  ring 
load  moving  with  a  constant  velocity  V,  as  shown  in  Fig.  1,  and  the  radial 
displacement  w  is  assumed  to  be  positive  inward.  The  axial  displacement 
and  the  pressure  and  axial  coordinates  are  denoted  by  TT,  p,  and  x,  respec¬ 
tively.  The  equations  of  motion  [5]  for  the  shell  are  then 


Eh  /a2u  v8w\  _ 

1  -  v ^  \3x^  a8x  /  8t2 


(1) 


and 


aEh' 


12(1  -  vc) 


a4w  ,  Eh  /w  8u\  —  a2w 

iFtr^'rvii;=pa-ph*«? 


(2) 


Next  the  following  variables 


2  2  — 

X  =  h/a  ,  a  =  E/p(l  -  v  )  ,  x  =  ax  , 


t  =  at/ a  ,  u  =  au  , 
p  =  pXE / ( 1  -  v2) 


w  =  aw 


(3) 


are  introduced  into  (1)  and  (2)  to  obtain  the  equations  of  motion  in  a 
dimensionless  form: 


82u  9w  32u 

57'  v^  =  ^ 


(4) 
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and 


8x4  8x  8t2 


For  the  moving  ring  load,  the  pressure  is  given  by 


p  =  pQ6(x  -  vt) 


(5) 


(6) 


where  pg  denotes  the  intensity  of  the  load,  6  the  Dirac  delta  function, 
and  v  the  nondimensional  speed  of  the  load.  The  actual  speed  V  of  the 
load  is  related  to  v  by 


V  =  av  .  (7) 

To  seek  the  steady-state  solution,  one  introduces  the  transformation 

z  =  x  -  vt  (8) 

and  sets  all  quantities  to  be  functions  of  the  new  variable  z.  This 
results  in 


(1  -  v2)u”  -  vw'  =  0  (9) 

.2  .  ? 

—  wlv  +  v  w"  +  w  -  vu'  =  Pq6(z)  (10) 


where  (  )'  =  d(  )/dz. 


o 


Operating  with  an  exponential  Fourier  transform  on  (9)  and  (10)  gives 


and 


where 


■>yP0 

6(1  -  v2)a 


(11) 


(12) 


f 


^!|4  -  V2!2  ♦  1  - 
12  b  * 


(13) 


(  I 


Here  £  is  the  transform  variable,  and  the  transform  of  a  quantity  (  )  is 
denoted  by  (  ),  Inverting  (11)  and  (12),  one  has 


and 


-  tVP°  r°  exp(-^z)  d£ 


=  [  exp(-t6z)d£ 

2tt  /  A 

J -00 


(14) 


(15) 


( 


! 
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Determination  of  Resonant  Speeds  for  Timoshenko  Equations 

Now  the  foregoing  integrals  are  well  defined  except  where  the 
velocity  is  such  that  A  has  a  double  root  on  the  real  axis.  The  condition 
for  this  to  be  so  is 


4 


v 


(16) 


For  these  values  of  v,  the  integrands  in  (14)  and  (15)  have  second-order 

poles  on  the  axis  of  integration,  and  for  such  a  situation,  the  integrals 

involved  do  not  exist  even  in  the  sense  of  a  Cauchy  principal  value.  It  will 

be  shown  later  in  the  analysis  that  such  nonexistence  of  a  solution  to  a 

physical  problem  implies  a  resonance  or  an  instability  in  the  sense  that  the 

displacements  everywhere  become  unbounded.  Values  of  v  for  which  (16)  is 

satisfied  are  the  resonant  speeds.  It  will  be  shown  that  there  is  one  root 

near  v  =  (X  /3)  '  and  two  near  v  =  1.  Actually,  there  are  other  roots  to 

(16)  also,  but  since  only  positive  real  roots  are  of  physical  significance, 

the  rest  are  discarded.  The  roots  can  be  tabulated  easily  as  functions  of 

v  and  X.  In  most  practical  problems,  X.  is  certainly  less  than  0.1  (i.  e.  ,  a 

thin  shell  is  assumed,  since  the  statement  X  <  0.  1  implies  that  the  thickness 

of  the  shell  is  less  than  0.  1  of  the  radius)  for  it  is  only  in  such  regimes  that 

the  shell  equations  are  valid.  For  values  of  X  >  0.  1,  it  is  doubtful  if  either 

Timoshenko's  or  Flvigge's  equations  are  valid.  Thus,  a  perturbation  solution 

in  X  is  a  valid  way  of  obtaining  the  approximate  roots  of  (16).  To  obtain  the 
2  1/4 

root  near  (X  /3)  ,  one  sets  v  =  v^j  to  be  given  by 


2 

VCR1 


=  «vQ  +*2vj  + 


(17) 
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where  c  =  X./n/3.  Substituting  in  (16)  and  equating  the  coefficients  of  the  like 
powers  of  c  ,  one  has 


VCR1 


2 

v  « 


4(1  -  v2)1/2 


+  • 


(18) 


Equation  (18)  is  valid  when  (1  -  v2)*^2  »  «  v2/4.  The  maximum  value 
of  v  is  0.5,  so  that  the  above  relation  is  always  true.  In  terms  of  physical 
velocity  V,  (18)  gives 


V2  Eh _ Eh2v2  +  -  .  . 

CR1  pa[3(l  -  v2)]1/2  6pa2(l  -  v2) 


(19) 


O 


The  first  term  of  (19)  is  the  square  of  the  bending  resonant  speed  given  in 
[4] ,  and  the  additional  terms  in  (19)  arise  from  the  consideration  of 
longitudinal  inertia  in  the  present  investigation.  It  is  apparent  from  (19) 
that  the  contribution  from  the  longitudinal  inertia  terms  is  indeed  small. 


To  obtain  the  roots  of  (16)  near  v  =  1.0,  one  sets  v  =  v^.^ 
given  by 


to  be 


vCR2  =  v0+<Zvl  +  -  •  •  <20) 

Substituting  from  (20)  into  (16)  yields 

2  ,  ,  2  2  .  .... 

vrn?  =  1  +  v  «  +  •  •  '  (21) 


( 


In  terms  of  physical  speed 


r  Zv2 
E  v  h 


CR2  7,  Z  ,  T7,  Z7 

p(l  -  V  )  3pa  (1  -  v  ) 


(22) 
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Another  resonant  speed  will  occur  when  the  constant  terms  in 
2  2 

&[i.e.,  1  -  v  /( 1  -  v  )]  defined  by  (13)  equal  zero.  This  situation  corre¬ 
sponds  to  a  double  root  at  £  *  0.  Setting  the  constant  term  to  zero  yields 
the  third  critical  speed 


(23) 


In  terms  of  physical  speed,  (23)  gives 


(24) 


2  2  2 

It  should  be  remarked  that  when  v  =  1,  i.e.  ,  V  =E/p(l-v),  there  is 
no  meaningful  solution  to  the  differential  equations  (9)  and  (10).  Thus, 
one  can  conclude  that  there  are  three  resonant  speeds,  given  by 


1. 


v2  _  _ Eh  _  EhZv2  . 

CR1  pa[3(  1  -  v2)]1/2  6pa2(l  -  v2) 


2. 


E 

P(1  -  v2) 


+  .  .  . 


3. 


E 

P 


The  first  resonance,  and  the  one  with  the  lowest  speed,  evidently  corresponds 
to  bending  resonance.  The  second  is  a  resonance  which  occurs  when  the  load 
speed  becomes  equal  to  the  plate-wave  speed  in  the  shell  and  is  an  axial 
vibration  resonance.  The  third  resonance  occurs  when  the  load  speed 
becomes  equal  to  the  bar-wave  speed  and  is  also  an  axial  vibration  resonance. 
It  is  not  clear  physically  as  to  why  the  bar -wave  speed  should  enter  at  all.  In 
fact,  for  load  speeds  of  the  order  of  magnitude  of  either  the  bar-  or  the 
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plate -wave  speed,  it  is  doubtful  that  the  equations  of  motion  without  the 
effects  of  shear  deformation  and  rotatory  inertia  are  valid. 

It  is  obvious  from  the  manner  in  which  the  instabilities  have  been 
obtained  here  that  these  critical  or  resonant  speeds  depend  only  on  the 
velocity  of  the  pressure  wave  and  not  upon  the  manner  in  which  the  pressure 
may  be  changing  with  z  above  and  behind  the  pressure  wave  front. 

Evaluation  of  Displacements  for  Timoshenko  Equations 

The  integrals  (14)  and  (16)  will  be  evaluated  next.  It  is  obvious  that 
the  integral  in  (14)  can  be  obtained  from  (15)  because 


u  =  -  - 


w' 


1  -  v4 


(25) 


Hence,  it  suffices  to  evaluate  (15),  which  in  explicit  terms  is 


w(z)  = 


lo 

2tt 


/oo 

00 


exp(-i£z) 


/ 4  -  v2|2  +  1  -  [v2/(l  -  v2)] 


di 


(26) 


First,  one  considers  the  case  when  v  <  v^^j.  The  roots  of  the 
denominator  of  the  integrand  are  given  by 


1 2  2v2  .2/,  v2  4\1/Z 

6  ) 


a  ±  ib 


(27) 


Since  v  <  v^j,  b  is  a  real  and  positive  number.  Then  the  roots  are 
located  at 


*  /l.  bl 

%  =  ±-=-±  1  — 

v*  ^ 


(28) 
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where 


and 


(29) 


(30) 


For  z  >  0,  it  suffices  to  integrate  along  the  real  axis  and  a  semicircular 
contour  in  the  lower -half  plane,  picking  up  the  contributions  from  the  poles 
at  |  =  ±a j/*'/*  -  ub y/*Jt  .  Performing  the  indicated  integration  results  in 


P0  •xp(-«b1/'/0  /  it, 

’ T5Tti .  „1,(1 .  v!„  1/it,  -  A(. -  A -  .Vi “'T*bl  T)  ‘  110 


(31) 


From  symmetry  considerations,  w  is  an  even  function  of  z;  hence 


Po  _ 

2^" [i  - 


v^/d  -  V 


exp(*bj/^)  /  **j  »» j  \ 

*» ■»(,  - ,»/(.  - 7i ■  A. ») w (‘» 101 T ~ T )  ’  ■£0 


(32) 


The  maximum  value  of  w  occurs  at  z  =  0  and  is  given  by 

w  p°  ln-iv2/i  ■vi)i1,z<v2),  l1/z 

max  2 [1  -  v2/(l  -  v2)]  1/2[1  -  v2/(l  -  v2)  -  v4/«  2]  1/2 

The  dynamic  amplification  factor  is  defined  here  as  the  ratio  of  the 
maximum  dynamic  deflection  to  the  maximum  static  deflection  (i.e.  ,  the 
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deflection  that  would  occur  when  v  =  0).  Hence,  the  dynamic  deflection 
(D.  A. )  is  given  by 


D.  A. 


r,  2,..  2.,  1/2,  2,  1/2,.  2.3/4 

[  1  -  V  /(I  -  V  )]  +  V  ft  (1  -  V  ) 

[l  -  V2/(l  -  V2)]  *^2[l  -  v2/(l  -  v2)  -  v^/c  2]  ^2 


(34) 


The  dynamic  amplification  evidently  becomes  unbounded  as  v  approaches 
one  of  the  critical  velocities. 

To  uncouple  the  effects  of  axial  motion  from  the  effects  of  radial 

motion,  it  suffices  to  neglect  the  term  v(9u/9x)  in  (5).  The  term  v(9w/9x) 

is  retained  in  (4)  so  that  the  axial  displacement  u  can  be  obtained.  The 

results  of  this  approximation  can  be  analyzed  from  the  theory  previously 

considered  in  the  present  paper.  The  neglect  of  v(9u/9x)  removes  the  term 
2  2-1 

-v  (1  -  v  )  in  the  equation  for  A,  viz.  ,  Eq.  (13).  Hence,  by  this  theory 
there  are  only  two  resonant  speeds  present,  and  they  are  given  by 


2 

VCR1  '  * 


V2  Eh 

CR1  a(l  -  v2)n/3 


(35) 


and 


2 

VCR2 


=  1 


or 


E 

P(1  -  v2) 


(36) 


The  first  resonant  speed  given  by  (35)  agrees  with  the  one  of  [4]  up  to  a 
2  1/2 

factor  of  (1  -  v  )  The  solution  for  the  displacement  now  becomes 


P0  exp(-zb  /_ 

2^(1  -(v</,2)]1/2V1 


COS 


z  >  0 


(37) 


and 
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z  <  0  (38) 


pQ  exp(zbjM)  /  zaj  _  zaA 

=  [I  -  "*  *  k‘ ,in  -5T/ 


where 


and 


(39) 


(40) 


and  v<vCRr 


The  maximum  deflection  again  occurs  at  z  =  0,  and  is 


w 


max 


,2*^  [l  -  (v2/*)]  1/2 


(41) 


The  dynamic  amplification  factor  is  obtained  from 


D.  A.  = 


(v*/.)]1'2 


(42) 


The  principal  differences  in  the  expression  given  above  are  in  the  neglect 
2  2 

of  the  v  /(I  -  v  )  terms.  The  first  resonant  speed  Vj-.Rj,  corresponding 
to  the  case  where  the  axial  coupling  is  not  neglected,  differs  from  the  case 


where  the  axial  coupling  is  neglected  by  terms  of  order  h/a.  But  because 
h/a  is  usually  quite  small,  this  approximation  is  certainly  a  valid  one.  Also, 
the  resonance  near  the  bar -wave  speed  disappears  if  the  axial  coupling  is 
neglected.  However,  since  the  shell  equations  probably  are  not  valid  in  this 
regime  of  the  velocities,  the  difference  is  not  of  much  importance. 


-14- 


Solution  Using  Flugge  Equations 


The  two  Timoshenko  solutions  will  be  compared  now  with  the  solution 
obtained  from  the  Flugge  equations  ([6],  pp.  209  -  219),  viz., 


(1  -  V  )u"  +  vw1  -  —  w'"  =  0 
4 


(43) 


2  2  . 

vu'  -  - —  u'"  +  - —  wlv  +  v  w 
4  4 


■+(1+t) 


w  =  p 


(44) 


Here  w  is  measured  positive  outward,  p  is  positive  when  it  is  a  suction 
rather  than  a  compression,  and  u  is  measured  positive  in  the  positive 
x-direction.  All  quantities  in  (43)  and  (44)  are  taken  to  be  functions  of 
z  =  x  -  vt,  and  the  notation  of  Flugge  is  changed  to  agree  with  that  of 
Timoshenko. 

Operating  with  Fourier  transforms  on  (43)  and  (44)  yields 


-|2(1  -  v2)u  -  +  L_  =  0 


(45) 


and 


-i^v  +L-42^r  +  -  v2|2  +^1  +^]w  = 


P0 


(46) 


where  p,  as  before,  has  been  taken  to  be  Pq6(z)  for  the  traveling  ring  load. 
Solving  for  u  and  w,  one  has 


i|[v  +  (*  2/4)$2]p0 

u  s  . . . — . ■—  ■ 

a*  d) 


(47) 
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and 


w 


iZd  -  v2)pq 

Zr(fj 


(48) 


where 


(1  -  v2) 


2*2 
v  £ 


(49) 


Again  it  suffices  to  study  only  w,  since  u  can  be  obtained  from  w  by  using 
(43).  Thus,  one  must  consider  the  integral 


w 


Pod  -O 


2ir 


r°°  2 

/  S  exp(-t|z)  . 

•'-CO 


(50) 


Determination  of  Resonant  Speeds  for  Fliigge  Equations 

As  before,  the  integral  will  not  exist  for  values  of  v  such  that  there 
are  second-order  poles  on  the  real  axis.  The  condition  for  this  to  be  so  is 


2.,  2.  „  v« ‘ 

r 11  -v  )+- 


(51) 


Equation  (bl)  is  a  higher  degree  equation  than  (16).  The  extra  root 
is  located  at  a  value  of  v  greater  than  1.0  and,  as  such,  is  not  of  much 
physical  importance  since  it  is  unlikely  that  the  shell  equations  without 
shear  deformation  and  rotatory  inertia  are  valid  for  values  of  v  in  this 
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o 


i 


range.  As  in  the  Timoshenko  equations,  there  is  a  root  nearv=  ‘J*  .  Setting 
v2  =  e  Xq  +  e  2Xj  gives  Xq  *  (1  .  v2)1/2  and  Xj  =  (-l/2)(v  +  v2).  Hence, 


2 

v 


-  |(v  +  V2)€  +  . 


(52) 


This  result  compares  well  with  the  result  from  the  use  of  the  Timoshenko 
equations  with  the  axial  coupling,  viz.  , 


2 

v 


(53) 


obtained  from  (18).  The  critical  physical  speed  is  given  by 


V2  = 


p[3(l 


E 

2TT 

v  )] 


h 

T7ZI 


1  (v+v2)E  h2 


Pd 


2,  2 

v  )  a 


(54) 


Again,  the  first  term  of  (54)  agrees  with  [4]  .  The  two  roots  greater  than 

2 

unity  are  very  close  together,  differing  from  unity  by  terms  of  order  c  and 

4 

from  each  other  by  terms  of  order  e  .  They  will  not  be  investigated  here, 

since  their  existence  is  doubtful,  because  when  the  equation  of  motion  was 

2 

derived,  terms  of  order  e  were  neglected.  Hence,  any  term  of  order 
2 

greater  thane  should  be  under  suspicion.  Dropping  terms  of  order  greater 
2 

thane  from  (51)  gives 


4..  2.2  2  2..  2.  2..  2...  2  2.  _ 

v(l-v)  +ve  v  (1  -  v  )  -  e  (1  -  v  )(1  -  v  -  v  )  =  0  (55) 

The  meaning  of  the  double  root  of  A' (£)  near  v  =  1  is  now  clear:  It  arises 

2 

from  the  factor  (1  -  v  )  common  to  all  terms  of  (55).  The  differential 
equations  themselves  predict  a  resonance  at  this  value  of  v  in  the  axial  mode. 
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2 

Canceling  the  factor  (1  -  v  )  yields 

v4(l  -  v2)  -  I  2[l  -  v2  -  (1  +  v)v2]  =  o 

The  corresponding  equation  from  the  Timoshenko  equations  from  (16)  is 

4..  2.  2,.  2  2.  n 

v  (1  -  v  )  -  *  (1  -  v  -  v  )  =  0 

The  root  near  v  =  ^*  is  unaffected,  but  the  root  near  v  =  1  now  becomes 


f 


(56) 


(57) 


2  .  ,  n  ,  ,2  2 

v  =  1  +  (1  +  v)  € 


(58) 


which  differs  from  the  results  of  the  Timoshenko  equations  by  a  term  of 
the  order  of  «  2.  The  Timoshenko  result  is  given  by  (21). 

Finally,  there  is  one  resonance  which  will  occur  if  the  constant  term 
in  (49)  disappears:  There  would  be  a  second-order  pole  at  the  origin  in  the 
integrand  for  w.  The  condition  that  the  constant  term  be  zero  is 


»  1 


(1  +€Z/4) 


2  ,  V2*  2 
v  +  — 


(59) 


This  value  of  v  will  always  be  slightly  less  than  unity. 

Evaluation  of  Displacements  for  Fliigge  Equations 

The  integral  for  w  will  now  be  evaluated.  The  procedure  and  results 
are  similar  to  those  for  the  Timoshenko  theory.  The  terms  of  order  greater 


O 
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2 

than  «  will  be  retained  until  the  end  of  the  analysis.  For  v  <  v^,  there 
are  four  simple  poles  located  at 


i  = 


(60) 


where 


2  [(1  -  €2/4p2)(l  -  v2/p2  +  «2/ 4)1  1/2  +  (v2/t  +  v«/202) 

*2  (l-£  2/402) 

,2  _  [q  -  «  2/4p2)(l  -  v2/p2  +  c  2/4)] 1/2  -  (v2/t  +  vt  /2p2) 
2  (l-*2/4p2) 


(61) 


(62) 


(  ' 


and 


(63) 


As  before,  the  integral  for  w  is  an  even  function  of  z  and  the  integral  for 
u  is  an  odd  function  of  z.  For  z  >  0,  the  use  of  a  semicircular  contour 
in  the  lower-half  plane  is  made  in  evaluating  (50).  The  result  of  the 
integration  is 


w 


exp(-zb£/^0 

2^  [(1  -  €  2/4p2Ml  -  v2/p2  +  «  2/4)  -  (v2/*  +  vt  / 202)2]  ^2 


Po 


x  (1  -«2/4p2)1/2 


(1  -  v2/p2  +  «  2/4) 


172 


a^  cos 


&)-.■“(?)]  ■ 


z  >  0 


•  (64) 


( 
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For  z  <  0,  one  substitutes  -z  for  z  in  (64).  The  maximum  deflection  again 
occurs  at  z  :  0  and  is  given  by 


2s/;  (i 


in .  (65) 

r"  «Z/4pZ)(l  -  vZ/pZ  +  «Z/4)  ,,*Z'Zi1/Z 


[(1 


(vZ/«  +v«/2pZ)Z]1/i 


The  maximum  static  deflection  is  obtained  by  putting  v  =  0  into  (65)  and 
is  given  by 


Po 


vno«  .tatic  =  (1  .  /  +  ,z/4)1/z  |[(l  -  «  z/4)(l  -  vz+«Z/4)]1/Z  -  v. /2 


J7T 


(66) 


The  dynamic  amplification  factor,  as  defined  earlier,  is 


p.  a.  = 


K1  -  <  2/*)u  -  v2 
2  / , - 2W i  Z/02  ,  2,,,!  iy 


1/2 


-  v«  /2 


1/2 


(1  -wi/pZ+«i/4)ll/z  |[(1  -«W)(1  -  vZ/pZ  +  «  2/4)]  Z,Z  -  <v*/«  +v«/2(J2)|i^ 


(67) 


Next,  the  powers  of  <  higher  than  first  will  be  neglected.  However, 

this  must  be  done  with  caution  since  the  order  of  magnitude  of  v  is  not  known 

precisely.  It  is  known  that  v  is  at  least  of  the  order  of  n/c  but  could  be  of 

higher  order  in  t  .  To  carry  out  this  procedure,  one  sets  v/n/c  =  0  and 

carries  0  in  all  terms,  assuming  0  is  of  order  one.  The  actual  procedure 

2 

then  is  to  carry  only  first  powers  of  £  ,  since  a  factor  of  c  is  already  taken 
out  by  setting  0  =  v/N It  .  The  results  are 


2  =  [(1  -  v2)  -  (2  -  v2)£02]1/2  +  02(1  -«02)  +  £  /2 


a2  = 


1  -  €0* 


(68) 
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(69) 


b2  -  Hi  -  y2>  -  <2  -  1/2  -  02d  -  <»2)  - 1  /a 

2-. 


.2x1/2 


w 


exp(-zb1M)(l  -  €0  )  P  2  2,  1/2 


Po_ _ 

2n/c  1(1  -  v2)  -  -  «[(2  +  v  -  v2)p2  -  20b] 


/  za2  za2\ 

XVa2COS^"+b2SinT"/  ’ 


z  >  0 


(70) 


w 


max 


D.  A.  = 


.  p0  |[(1  .  v2)  -  .(2  .  v2)02)1/2  +  »2  *  [(v/2)  -  02].|1/2(1  -<*2) 

■  *E  <77  ,rr.A"* |„  :  V2,  -  /.  .[,2  ♦  v  ♦  vV  -  2,‘]  |*»  '  ’ 

(1  -  v2)(l  -  «»2)[<1  -  v2)1/2  -  v./2]1/2 

o  -v2-  ./,■'>  -,»T7* 02-[(v/2,.02i.|'/2  (72> 

Response  in  the  Supercritical  Range 


For  values  of  v  between  the  first  and  second  resonances,  given  by 
(19)  and  (22),  respectively,  the  integral  (15)  defining  w  can  also  be  evaluated. 
In  this  case,  there  will  be  four  simple  poles,  all  on  the  real  axis  and  located 
at  |  =  ±Qj,  ±02'  °nlY  tl'e  Timoshenko  solution  will  be  discussed  for  this 
case.  The  poles  are  given  by 


and 


v2  + 


1/2 


(73) 


(74) 
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The  integral  (15)  then  is  defined  in  the  sense  of  a  Cauchy  principal  value 
and  is 


p0  1 

'sin  (agz) 

8  [v4-«*(i  -v*/e*)]1/z 

a2 

sin  (eTj*) 
~  , 


z  >  0 


(75) 


where  *  (2 /«  )aj  and  ag  =  (2/*  Jag-  The  solution  for  z  <  0  is  obtained  by 
setting  -z  for  z  in  (75). 

The  physical  meaning  of  the  solution  (75)  is  not  at  all  clear.  The 
solution  predicts  that  both  the  deflection  and  the  slope  are  zero  under  the 
load  and,  thus,  supplies  an  utterly  improbable  result;  such  a  result  does 
not  appear  to  be  physically  realizable.  It  is  felt  that  for  values  of  v  in 
this  range,  an  initial  value  problem  should  be  solved,  since  it  is  known 
[7]  that  the  steady-state  problem  is  an  artificial  one.  Values  of  v  should 
really  be  derived  from  a  transient  solution  by  allowing  time  to  become 
unboundedly  large.  This  problem  is  beyond  the  scope  of  the  present 
investigation.  For  a  related  analysis,  the  reader  is  referred  to  [8] . 
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SECTION  IV 
SUMMARY 

The  stability  of  a  circular  cylindrical  shell  subjected  to  a  moving  ring 
load  with  a  constant  velocity  has  been  examined  in  detail  when  both  longitu¬ 
dinal  and  transverse  inertia  effects  are  included,  using  both  the  Timoshenko 
and  the  Flligge  equations.  It  has  been  found  that  the  resonance  due  to  bending 
occurs  at  a  lower  velocity  than  that  reported  in  [  1  ] .  The  resonance  due  to 
bending  is  concealed  in  the  solution  of  [  l]  because  of  the  assumed  presence 
of  viscous  damping  in  the  radial  direction.  However,  a  bending  resonance 
occurs  at  a  much  lower  velocity  than  that  in  [  l]  and  is  of  importance  for 
design  of  aerospace  vehicles.  The  instabilities  have  been  obtained  in  the 
present  paper  using  Fourier  transforms  without  obtaining  the  entire  solution. 
The  response  of  the  shell  is  obtained  in  the  subcritical  range  for  Timoshenko 
and  Fliigge  equations.  Expressions  for  the  dynamic  amplification  in  the 
deflection  have  been  given.  The  critical  speeds  obtained  by  using  Flligge 
equations  differ  by  only  a  small  amount  from  those  given  by  the  Timoshenko 
equations. 

A  formal  solution  for  the  response  in  the  supercritical  range  has  been 
obtained  for  the  Timoshenko  equations.  The  formal  solution  does  not  appear 
to  be  physically  meaningful.  Therefore,  it  is  felt  that  an  initial  value  problem 
should  be  solved,  including  the  effects  of  shear  deformation  and  rotatory 
inertia,  for  values  of  velocity  in  the  supercritical  range. 


I  ( 
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